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The generation of an early kination dominated (KD) era within a tracking quintessential model is investi- 
gated, the relic density of the Weakly Interacting Massive Particles (WIMPs) is calculated and we show that 
it can be enhanced with respect to (w.r.t) its value in the Standard Cosmology (SC). By adjusting the param- 
eters of the quintessential scenario, the cold dark matter abundance in the universe can become compatible 
with large values for the annihilation cross section times the velocity of the WIMPs. Using these values and 
assuming that the WIMPs annihilate predominantly to p + p~, we calculate the induced fluxes of e cosmic 
rays and fit the current PAMELA and Fermi-LAT data. We achieve rather good fits in conjunction with a 
marginal fulfillment of the restriction arisen from the Cosmic Microwave Background (CMB). 



1 Introduction 

A plethora of recent data | l j indicates that the two major components of the universe are the Cold Dark 
Matter (CDM) and Dark Energy (DE). The DE component can be explained with the introduction of a 
slowly evolving today scalar field, q, called quintessence whereas WIMPs, x, are the most natural candi- 
dates to account for CDM. In this talk, which is based on Ref. Q, we reconsider (Sec. [2J the creation of 
an early era dominated by the kinetic energy of q in the context of a tracking quintessential model 0. We 
show (Sec. [3} that if \ decouples from the cosmic fluid during this era, its relic abundance, tt x h 2 can be 
significantly enhanced w.r.t. its value in SC. This enhancement of tt x h 2 assists us to interpret, through 
the x's annihilation in the galaxy (Sec. |4), the reported 015] excess on the positron (e + ) and/or electron 
(e — ) cosmic-ray (CR) flux under the assumption that x's annihilate predominantly into p + pT (Sec.|5j. 
Throughout, the subscript [I] is referred to present-day values [to values at the onset of our scenario] and 
pi = Pi/p c o (i = g,R and M) where p c0 = 8.1 • l(T 47 /i 2 GeV 4 with h = 0.72. 

2 The Tracking Quintessential Model 

The quintessence field, q, of our quintessential scenario (QS) satisfies the equation: 



q+ZHq+dV/dq = 0, with V = M i+a /q a + b H 2 q 2 /2 and H/H = H= ^p q + pR + pM- (1) 

Here H is the Hubble parameter, dot denotes derivative w.r.t the cosmic time t, p q = q 2 /2 + V, pr ~ 
Pro exp ( — 4x) and pm = Pmo exp (— 3x) is the q, radiation and matter energy density respectively, x = 
ln(R/Ro) is the logarithmic time and R is the scale factor of the universe. 
We impose on our QS the following constraints: 

(a) Initial Domination of Kination. We focus our attention in the range of parameters with Vl q \ = 
^q(^i) = 1 where il q ~ p q /{p q + Pr + Pm) is the quintessential energy-density parameter. 

(b) Inflationary Constraint. Assuming that the power spectrum of the curvature perturbations is generated 
by an early inflationary stage, we impose the bound Hi < 1.72 • 10 56 . 

(c) Nucleosynthesis (BBN) Constraint. At the onset of BBN, Xbbn = —22.3, p q is to be sufficiently 
suppressed compared to pr, i.e., J6) ^(xbbn) < 0.21 at 95% confidence level (c.l.). 
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FIG. 1 The evolution of log pi with i — q (gray [dark 
gray] line), R+M (light gray line) and A (thick line) 
as a function of x for qi = 0.01, a = 0.5, b = 0.2, 
Ti = 10 9 GeV, M = 4.8 eV and Hi = 4.7 ■ 10 51 
[Hi = 2 • 10 53 j - note that pn +M = PR + pm- We ob- 
serve that although the used Hi 's differ by two orders 
of magnitude, both p q 's reach p\ highlighting thereby 
the insensitivity of our QS to the initial conditions. 
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(d) DE and Cosmic Coincidence Constraint. These two requirements can be addressed if we demand Q 
tt q o = p q o = 0.74 and d 2 V(i; = 0)/dq 2 ~ H%. 

(e) Acceleration Constraint. Successful quintessence has to account for the present-day acceleration of 
the universe, i.e. fl] -1.12 < w q (0) < -0.86 (95% c.l.), where w q = (q 2 /2 - V)/(q 2 /2 + V) is the 
barotropic index of the q-field. 

Solving Eq. (fl} we find that q undergoes four phases during its evolution - see Fig. Q] For x < Xkr 
the universe and q is dominated by q/2 3> V and therefore we get a KD era, during which q is set in 
anharmonic oscillations for b > 0. In particular, q develops extrema at 



For x > Xkr the universe becomes initially radiation and then matter dominated whereas p q is dominated 
initially by q/2 and then by V. As x approaches the system in Eq. (fl} admits [3 1 a tracking solution since 
the energy density of the attractor: 



tracks pM until it outstrips and dominates the current expansion of the universe. It can be shown f 2 ,.3 ] that 
b > ensures the coexistence of an early KD phase with the achievement of the tracking solution in time. 
Moreover, for a < 0.6 the requirement 2-(e) can be marginally fulfilled, too. 

3 The WIMP Relic Density 

The relic density of a WIMP x> ^ x /i 2 , with mass m x is calculated using the formula: 

Q x h 2 = 2.748 • 10 8 (n x o/s ) (m x /GeV) , where h x + 3Hn x + (av) {n 2 , - rtf) = (4) 

is the Boltzmann equation which governs the evolution of the %'s number density, n x . Also, n° q oc 
x 3/2 g-i/a w j tn x _ T/m x is the equilibrium configuration of n x , (av) is the thermal-averaged cross 
section times the velocity of %'s and s oc T 3 is the entropy density with T the temperature. 

The decoupling of x from the cosmic fluid during the QS and SC is visualized in Fig. 2-(a). We observe 
that, in both cases, the current p x /s follows p x q /s and at some x = Xp, p x /s dominates over p x q /s and 
remains constant until today. For the selected m x and (av) we obtain an enhancement of il x h 2 within QS 
w.r.t its value in the SC, Q x h 2 \ sc , since Q x h 2 = 0.11 whereas Sl x h 2 \ sc = 0.0045. This enhancement 
can be further analyzed, by defining AQ X = Q x h 2 / Q x h 2 \ sc — 1. The behavior of Afi x as a function 
of the free parameters of the QS can be inferred from Fig. 2-(b). For b = we obtain a pure KD era and 
Af2 x increases when m x increases or (av) decreases. For 6^0, Afl x depends crucially on the hierarchy 
between Xp and x cx t- As m x increases above 0.1 TeV, Xp decreases and moves closer to x cx t and Afl x 
decreases with its minimum Afi x | • occurring at m x = m x \ . which corresponds to Xp m ~ T ex t- 




(2) 




(3) 
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FIG. 2 (a) The evolution as a function of': of the quantities p x q /s = m x n x q /s (dotted line) and p x /s — 



(thick [thin] solid lines) for the QS [SC]; (b) Af2 x versus m- 



(av) = 10~ 6 GeV" 2 [(av) 
and b = 0.32 (dotted lines). 



for the QS with a = 0.5, Hi = 6.3 • 10 53 , Ti = 10 9 GeV, 



10 GeV ] (gray [light gray] lines) and b = (solid lines), b = 0.15 (dashed lines) 



4 e ± -CRs From WIMP Annihilation 

Residual x's annihilations in the galaxy induce a e + flux per energy at Earth which is given by 

v e + is the velocity of e + , p Q = 0.3 GeV/cm 3 is the local CDM density, b(E) = E 2 /(GeVt E ) with 
ts = 10 16 s is the energy loss rate function and dN e + /dE e + denotes the energy distribution of e + 's per 
X annihilation and can be found in Ref. (7). Also, /(Ad) is the dimensionless halo function which fully 
encodes the galactic astrophysics and can be read off from Ref. JU 

There are three sources of uncertainty in our computation: the CDM distribution, the propagation of \ 
annihilation products and the astrophysical backgrounds. We adopt (a) the isothermal halo profile, to avoid 
troubles with observations on 7-CRs; (b) the MED propagation model, which provides the best fits to the 
combinations of the two data-sets; and (c) commonly [8| assumed background e fluxes normalized with 
the Fermi-LAT data. Adding the latter contributions to the one in Eq. (f5J) we get the total fluxes, $ e ± . 

In order to qualify our fittings to the experimental data, we define the % 2 variables as follows: 

( Fm-*m? „,,, „ _ r$ e+ /($ e+ +$ e -) andiV A = 7 for PAMELA, 



2 \ - y Ai ~ J -Ai) - . F = V/l* e +t* e -J mm iM A = 1 101 rAivinL.A, 
*>A (AF^ bs ) 2 ' A \ J B; 3 +($ e ++$ e -) and Na = 26 for Fermi LAT, 

where i runs over the data points of each experiment A, "obs" ["th"] stands for measured [theoretically 
predicted] values. The best fits to the combined experimental data can be achieved with m x ~ 1.28 TeV 
and (av) ~ 1.95 ■ 10~ 6 GeV" 2 resulting to (% 2 + xj)/d.o.f = 24/31. 

5 Results 

To systematize our approach, we can define regions in the m x — (av) plane which are favored at 95% c.l. 
[99% c.l.] by the various experimental data on the e ± -CRs demanding 



X 2 <X 2 L„+6 X 2 < X 2 1 +9.2 



with % 2 



X{ for PAMELA, 

Xi + X\ for PAMELA and Fermi LAT, 



where % 2 1 . can be extracted numerically by minimization of x 2 w.r.t m x and (av). 

The large (av)'s which are required in order to fit the experimental data on e ± -CRs are to be consistent 
with a number of requirements so as the interpretation of the data on e ± -CRs via CDM annihilation in the 
galaxy is fully acceptable. All in all, we impose the following constraints: 
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FIG. 3 Restrictions in the m x — (av) plane for a — 0.5, T\ = 10 9 GeV and several b's and Hi's indicated in 
the graph. The light gray shaded areas are allowed the constraint 5-(a), the sparse black hatched area is preferred at 
95% c.J. by the PAMELA data and the dense black [red] hatched areas are preferred at 95% c.l. [99% c.J.j by the 
PAMELA and Fermi-LAT data. Regions above the black solid, dashed, dot-dashed and dotted lines are ruled out by 
the upper bounds on (av) from the constraints 5-(b), 5-(c) 5-(d) and 5-(e), respectively. The best-ht (m x , (av)) for the 
combination of PAMELA and Fermi-LAT data which is consistent with all the constraints is given in the table. Shown 



are also Q, x h 



and several b's and Hi 's (and the resulting Tkr 's) leading to Q, x h ~ 0.11 in our QS. 



(a) Constraint from the CDM abundance £TJ: 0.097 < Q, x h 2 < 0.12. 

(b) BBN constraint (SI: (av) < 8.6 • 10~ 5 GeV" 2 (m x /TeV) . 

(c) CMB constraint fifflj: (av) < 1.3 • 10~ 6 GeV~ 2 (m x /TeV) . 

(d) Constraint from the j-CRs ^771/ : (av) < 4 • 10~ 6 GeV~ 2 for the isothermal halo profile. 

(e) Unitarity constraint: (av) < 8ir GeV~ 2 (m x /GeV) 2 . 

Imposing all the constraints above we can delineate our findings in the m x — (av) plane as in Fig. [3] 
A simultaneous interpretation of the e ± -CR anomalies consistently with the various constraints can be 
achieved in the regions where the gray shaded areas overlap the lined ones below the dashed lines. We 
observe that part of the region favored at 99% c.l. by PAMELA and Fermi LAT is allowed. The best-fit 
(m x , (av)) - with x 2 /d.o.f = 33/31 - which saturates the most stringent (CMB) bound is arranged in the 
Table of Fig. [3] We remark that the requirement 5-(a) is violated within SC but can be met by adjusting the 
parameters of the QKS. In all cases we obtain Tkr < 0.02 GeV with Tkr being the transition temperature 
to the conventional RD era. It remains the construction of a particle model with the appropriate couplings so 
that x's annihilate into fi + fi~ with the desired (av) derived self-consistently with the (s)particle spectrum. 
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